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Abstract 

We study the modified theory of gravity in Friedmann Robertson 
Walker universe composed of several perfect fluids. We consider the power 
law inflation and determine the equation of state parameters in terms of 
the parameters of modified gravity's Lagrangian in the early universe. We 
also discuss briefly the gravitational baryogenesis in this model. 

1 Introduction 

Modified theory of gravity, constructed by adding geometrical correction terms 
to the usual Einstein Hilbert Lagrangian has been used to study the inflationary 
epoch and the subsequent reheating stage in the early universe [TJ [2] . In this 
view the effective Lagrangian, L, in the early universe includes higher order 
curvature terms and inflation may be a natural result of this theory. In pQ, 
it was shown that involving a term proportional to the square of the scalar 
curvature, i.e., L = R + aR 2 , results in a quasi-de Sitter expansion. In this 
model the Hubble parameter decreases slowly for large a before going into an 
oscillation phase which can reheat the universe [5]. In principle, one can assume 
that the effective Lagrangian is a function of the scalar curvature L — f(R) [3]. 
The time dependence of the scale factor depends on f(R), e.g., if one chooses 
f(R) = R + aR 2 + f3R 3 , instead of an exponential like inflation, he may obtain 
a power law expansion for the universe Some of particle physics problems 
such as the hierarchy problem [5] and baryogenesis [H] can also be studied in 
the framework of modified theory of gravity. Recently, the modified theory of 
gravity which is able to describe the present cosmic acceleration [7] without 
involving exotic dark energy [5] has attracted more attention. 

In this paper we aim to study the universe in an era when the expansion can 
be described by the scale factor a(t) oc t x . This assumption is consistent with 
most parts of the (Friedmann Robertson Walker expanding) universe history, 
when one of the components such as (dark) matter or radiation dominates. But 
in the framework of the usual theory of gravity, and in the presence of ordinary 
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matter and radiation, the power law expansion can not describe the inflation- 
ary and accelerating expansion phases. Besides, the gravitational baryogene- 
sis proposal fails to determine the baryon asymmetry in radiation dominated 
era. In this paper we consider the modified theory of gravity, characterized 
by f(R) = '^2 li AiR n \ in the early universe. We assume that the universe is 
composed of different perfect fluid components with constant equation of state 
parameters. Physical parameters of the system such as radiation temperature 
and the equation of state parameters of the fluids are determined in terms of 
rij's and A^s. At the end we use our results to determine the gravitational 
baryogenesis in a universe composed of two components in the framework of 
modified theory of gravity. Through the paper we use the units c = fee = !• 



2 Power law expansion in modified theory of 
gravity 

We consider the modified theory of gravity described by the action 

S = J {^f{R) + L m )^—gd i x, (1) 

where L m is the Lagrangian corresponding to the matter such as radiation, 
baryonic matter, dark matter and so on. The geometry of the universe which is 
assumed to be spatially flat, homogeneous and isotropic is described by Fried- 
mann Robertson Walker (FRW) metric 



ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ). (2) 

jft) 
i(t)- 



a(t) is the scale factor in terms of which the Hubble parameter is H — The 



Ricci scalar is 

R = Q(H + 2H 2 ). (3) 

By variation of the action with respect to the metric, we obtain 

g^af(R) - V^ v f'(R) + f'R pv - = SttGT™, (4) 

where T™ v are the energy momentum tensor components of matter fields which 
behave like a perfect fluid. The energy density, p, and the pressure, P, may be 
derived from |4|: 

^Gp = &-+^f\R)R-£f{R) 
la a 



87tGP = -i^ + {2 ^ + -)f(R)-2-f(R)R-r(R)R 2 - 

f"(R)R. (5) 

The energy conservation relation, 

p + 3H(P + p) = 0, (6) 

is not independent of (0, but is required for consistency. We assume that the 
perfect fluid is effectively composed of non-interacting components, although 
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each component may include subcomponents whose interactions are rapid com- 
pared to the expansion rate keeping them in thermal equilibrium. In this way 
the total energy density and the total pressure are given by 

i i 

Each component satisfies 

Pi + 3H(Pi + Pi ) = 0. (8) 

For time independent equation of state parameters (EOS) (denoted by 7i's), the 
solution of J5|), in terms of the scale factor, is 

Note that even for constant 7^'s, the EOS parameter of the universe, 

7=- = ^^, (10) 
9 2Z l 9 l 

is time dependent. To study (0, let us take 

f(R) = Y,fi(R), (ii) 

i 

such that fi satisfies 

2 a a 

2 a z a a 

-fi'(R)R. (12) 

In this way if (fl2|) is satisfied, then ([5|) will be also satisfied. In this paper we 
restrict ourselves to the models characterized by [5] 

f(R)=J2AiR ni . (13) 

i 

AiS and n^'s are real constants (n^s are not restricted to integer numbers). 
Note that f(R) broken power law models such as [TU] 

f(R) = R-m 2 1{ ^> 



+ 1' 

where to, n > 0, ci, ci are real numbers, can be casted to (fT3"f at high curvature 
limit R >• m 2 . 

f(R) in |T|) must satisfy some stability conditions [H]. For the model ((T|), 
the condition d J^p > 0, which in our model reduces to 

J2M^-±)A l R ,H ~ 2 >o, (14) 
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is necessary for classical stability of FRW solution in high-curvature regime. 
In the context of quantum mechanics this condition ensures the absence of 
tachyonic scalarons [T]. We also require that ji? > or 

J^niAiR 71 *- 1 >0, (15) 

i 

to prevent the graviton from turning into a ghost |llj . 

In the following We use the ansatz a(t) oc t x for the scale factor. This ansatz 
is allowed when the number of (dominant) fluid components is the same as the 
number of the terms in f(R). In general, A depends on EOS parameters of the 
components of the perfect fluid, 7i's, and rij's. Now let us find the conditions 
required for consistency of a(t) oc t x with l|12p. By substituting the Ricci scalar, 
R = 6 A and © into $T2§ we obtain 

/ \ / t \ -3A(7i + l) 

(2A-l)^- 1 (6A)^T*(3-A)-2n?-l+2AjA i t- 2n < = l&irG/n(jt )^—J 

(16) 

The above equation is true for all t, provided that 

3A m rii 

V«,J, (17) 



2 l+ 7l 1 + 7j 

and 



Ai = PSrt* — (18) 

<5 4 



where 



ft = (n,(-A + 3) - 2n 2 + 2A - l) (6(2A - l))" i_1 A™\ (19) 

Hence if the EOS parameters of the fluid components satisfy (fTT)) . and Aj's 
are given by (fT%|) . then in /(i?) gravity specified by (fl"3j) . a(t) oc t x can be 
considered as the scale factor. If m — 0, then we must take 7; = —1 describing 
a cosmological constant corresponding to the vacuum energy. For positive A, 
n j < leads to 7^ < — 1 describing a phantom like dark energy component. 

The above discussion can be generalized to the case that A is a slowly varying 
function of time, A = A(t), X(t)t <C X(t) [S]. Time independence of Ai requires 

Pi R Qi 



where Q t is defined through Q t = R 2 - 671* [H + H 2 j R + 6m(m - 1)HR. By 
7i = —1 — 3 £ p and (|20j) . we can determine the time dependent EOS parameters 

3X-2n t 4(4A-l)n 2 + (4A 2 -17A + 2)n l + llA-8A 2 Ai 
'" ~ ~ + 3A(2A - l)(2n 2 + (A - 3)m + (1 - 2A)) T 

(21) 
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Up to O fx)' ^ ne re l a tion between ji's may be obtained in a compact form 
7i(*) = ~7i(*) + (--l), Vi,j 

= - 1+ mr < 22 > 

It is worth to note that in a one component fluid, with EOS parameter 7, we 
have 7 = |j — 1, hence if the universe is approximately filled with radiation (i.e 
7=|) then A = § and i? = 3n ( s p^) ■ Therefore for n 7^ 1 we have i? ^ and 
i? 7^ while in Einstein theory of gravity we obtain R = 0. This is the note 
used in [5] to show that in modified gravity the gravitational baryogenesis may 
occur even in radiation dominated epoch. 

Now assume that one of the thermal fluid components (e.g. the radiation 
component) has temperature T, i.e., the subcomponents of this component have 
nearly common temperature T. We assume also that the density of each com- 
ponent is given by 

Pi = eiT e \ (23) 

This is allowed when T is proportional to a power of the scale factor. It is 
worth to note that T may not be the temperature of other noninteracting fluid 
components. and (|18jl can be used to obtain the time dependence of T 

T=(sAyt^. (24) 
This equation holds for each i, therefore 



m n 



_ 2 

Besides, for Vi,j we must have 



(25) 



^AA =\^A j \ . (26) 



To elucidate our results, as an example, we assume that the universe is approx- 
imately composed of a radiation component (denoted by the subscript TV) and 
a non thermal component with EOS parameter uj [12 ]. These non-interacting 
components satisfy the energy conservation relation 

p n + AHp n = 
p u + 3H(l+u)p u = 0. (27) 

The time derivative of the ratio of these components, denoted by r := is 

r = 3Hr (u> - ~) (28) 

if lj > ^, r > and p u component decreases more rapidly than radiation 
component and the universe will become radiation dominated. Following our 
previous discussions let us take 

f(R)=A n R n * + A u R n ». (29) 
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Stability conditions require that the parameters of the model satisfy 



n n (nn - l)A n R n ^ + n w (n w - > 0, (30) 

and 

rmAnR^- 1 + n^R^- 1 > 0. (31) 
If the scale factor is given by a(t) oc t x , then 

n n = 2\, n w = ^(l+w)A. (32) 

As a result, in this model, uj can be expressed in terms of Uj's: lo = — 1. 
The stability conditions is satisfied by choosing appropriate parameters for the 
model, e.g., if we take n u = 1, which results in ui — A — 1, (|29|) reduces to 
f(R) = A U R + A-jiR 1171 (the stability of models including this specific case was 
discussed in [13j ) and the stabilities conditions become 

n n (nn ~ l)A n R n ^ > 

nnAnR^- 1 + A u > 0. (33) 

For positive curvature and for A's belonging to the domain 0.5 < A < 0.64, An 
becomes a positive real number (see (fT!?)) and Besides, 5^ — implies 

that A^ is positive. Therefore the stability conditions are satisfied. Note that in 
[H] it was proposed that f(R) models with B < 0, where B is defined through 



d 2 f 



dR fdluH 



B= ^i^{i^) > (34) 

dR v ' 

are unstable to linear perturbations at high curvature, this leads us to take 
A > 0.5 which is in in agreement with 0.5 < A < 0.64 proposed above. 
Following and (TT%|) . the energy densities can be obtained as 

Pi = M;i- 3A(1+7i) . (35) 
For radiation component this yields 

Pn = ^(-10A 2 + 8A - 1)A 2A (6(2A - l)) 2 ^ 1 A n t^\ (36) 

07TG 

On the other hand the temperature of the radiation component is given [15j 

Pn = enT\ (37) 

2 

where e-ji — f^ff* and g+ is the total degrees of freedom of effective massless 
particles contributing in radiation component. Hence, in (|23|) . 8-jz — 4. By 
using (|37[) and ([33)1 we find out 

T= fS*£*y t -*¥>. (38) 

By substituting (|38|) into ([35| one gets 



T^ "^ 4 - (39) 



6 



By considering the assumption (|2"5|) we find 



e u = 5 U A U ( -p- ) " n , U = 4^. (40) 

If at a time denoted by t = tuo , the energy density of radiation component 
becomes equal to the other component: 

Pu(tRD) = Piz(t RD ) = e-RT^Di ( 41 ) 

where Tro is the radiation temperature at tRo, then the solution of |2T|) may 
be written as 

P» - e n T RD (J!L\ . (42) 

The temperature Tro can be determined in terms of the parameters of the 
model: 



Trd = e K 4 {SMMn n - nu ) (5 n A n y^-m) . (43) 

As an application, we can use these results to study a simple gravitational 
baryogenesis model [T^] , in the context of modified theory of gravity in the early 
universe composed approximately of radiation and a nonthermal component 
described by |2T|) . The key ingredient in this theory is the coupling of derivative 
of the Ricci scalar curvature and the baryon number current 



£ 

A 2 " 



J d 4 x^(d^R)J^ (44) 



where A is a cutoff characterizing the scale of the energy in the effective theory 
and s = ±1. (144|) dynamically violates CPT giving rise to the baryon asymme- 
try. In a universe with spatially constant R, to obtain the chemical potential 
for baryon (pb) and antibaryons {pg), we use 

±^R)J» = ^R(n B -n s ), (45) 



where tib and n§ are the baryon and antibaryon number densities respectively. 

2sR 

A 2 



Therefore there is an energy shift, for a baryon with respect to an an- 



tibaryon. We can assign a chemical potential to baryons: /ib = — Pb = — 
So, in thermal equilibrium there will be a nonzero baryon number density given 
by: 

n b = n B - ng = \ tx — + {—) A j , (46) 

where g b ~ 0{\) is the number of internal degrees of freedom of baryons. The 
entropy density of the universe is given by s — ^-g s T 3 , where g s ~ 106 indicates 
the total degrees of freedom for relativistic particles contributing to the entropy 
of the universe 15J. In the expanding universe the baryon number violation 
decouples at a temperature denoted by To and a net baryon asymmetry remains. 
The ratio ^ in the limit To S> m b (nib indicates the baryon mass), and To ^> p b 
is then: 

^~_ £ 2^A| T (47) 
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Introduction of £ gives us the possibility to choose the appropriate sign for rib- 
In our model R is determined as 

R = -6n K (l - n n ) (j^A . (48) 

To derive (|48|) we have used the fact that the decoupling time is 

By putting (|48p into (|47|) we can determine the baryon asymmetry 
n B 0.02 

Note that (f5U|) is valid for radiation dominated epoch as well as for non-radiation 
dominated era provided that the EOS parameter of non radiation component 
satisfies (|32|) . 
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